Ordering of the three-dimensional Heisenberg spin glass in magnetic fields 
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Spin and chirality orderings of the three-dimensional Heisenberg spin glass are studied under 
magnetic fields in light of the recently developed spin-chirality decoupling-recoupling scenario. It is 
found by Monte Carlo simulations that the chiral-glass transition and the chiral-glass ordered state, 
which are essentially of the same character as their zero-field counterparts, occur under magnetic 
fields. Implication to experimental phase diagram is discussed. 



In the studies of spin glasses, much effort has been 
devoted to the properties under magnetic fields. Unfor- 
tunately, our understanding of them still has remained 
unsatisfactory fij]. Most of the numerical studies have 
focused on the properties of the simple Ising model, es- 
pecially the three-dimensional (3D) Edwards-Anderson 
(EA) model. While the existence of a true thermody- 
namic spin-glass (SG) transition has been established for 
this model in zero field, the question of its existence or 
nonexistence in magnetic fields has remained unsettled. 

If one tries to understand real experimental SG order- 
ing, one has to remember that many of real SG materials 
are more or less Heisenberg-like rather than Ising, in the 
sense that the random magnetic anisotropy is consider- 
ably weaker than the isotropic exchange interaction [|],|| . 
Numerical simulations have indicated that the isotropic 
3D Heisenberg SG with finite-range interaction does not 
exhibit the conventional SG order at finite temperature 
in zero field jj]-^). Since applied fields generally tend 
to suppress the SG ordering, a true thermodynamic SG 
transition is even more unlikely under magnetic fields in 
case of the Heisenberg model. 

Experimentally, however, a rather sharp transition-like 
behavior has been observed under magnetic fields in typ- 
ical Heisenberg-like SG magnets, e.g., canonical SG like 
AuFe and CuMn, although it is not completely clear 
whether the observed anomaly corresponds to a true ther- 
modynamic transition The situation is in con- 
trast to the zero-field case where the existence of a true 
thermodynamic SG transition has been established ex- 
perimentally [Q. Set aside the question of the strict 
nature of the SG "transition" , it is experimentally ob- 
served that a weak applied field lowers the zero-field 
SG transition temperature rather quickly For 
higher fields, the SG "transition" becomes much more ro- 
bust to fields, where the "transition temperature" shows 
much less field dependence ji]]6jQ] . Such behaviors of the 
SG transition temperature under magnetic fields T g {H) 
were often interpreted in terms of the mean-field model 
Indeed, the mean- field Sherrington-Kirkpatrick 
(SK) model with an infinite-range Heisenberg exchange 
interaction with weak random magnetic anisotropy ex- 



hibits a transition line similar to the experimental one 
||: i.e., the so-called de Almeida-Thouless (AT) line 
H (x (T g (0) — Tg(H)) 3 / 2 in weak-field regime where the 
anisotropy is important, and the Gabay- Toulouse (GT) 
line H cx (T s (0) — Tg(H)) 1 ^ 2 in strong-field regime where 
the anisotropy is unimportant. Nevertheless, if one notes 
that the true finite-temperature transition under mag- 
netic fields, though possible in the infinite-range SK 
model, is unlikely to occur in a more realistic finite-range 
Heisenberg model, an apparent success of the mean- 
field model in explaining the experimental phase diagram 
should be taken with strong reservation. 

Recently, one of the present authors has proposed 
a scenario, the spin-chirality decoupling-recoupling sce- 
nario, aimed at explaining some of the puzzles concerning 
the experimentally observed SG transition in zero field 
In this scenario, chirality, which is a multispin vari- 
able representing the sense or the handedness of local 
noncoplanar spin structures induced by spin frustration, 
plays an essential role. In a fully isotropic Heisenberg 
SG, in particular, this scenario claims the occurrence of 
a novel chiral-glass ordered state in which only the chi- 
rality exhibits a glassy long-range order (LRO) while the 
spin remains paramagnetic. At the chiral-glass transi- 
tion, among the global symmetries of the Hamiltonian, 
only the Z2 spin reflection (inversion) symmetry is bro- 
ken spontaneously with keeping the SO(3) spin rotation 
symmetry preserved. Note that this picture entails the 
spin-chirality (or SO (3) — Z 2 ) decoupling on long length 
and time scales: Namely, although the chirality is not 
independent of the spin on microscopic length scale, it 
eventually exhibits a long-distance behavior entirely dif- 
ferent from the spin. Such a chiral-glass transition was 
indeed observed in zero field in a recent Monte Carlo 
simulation by Hukushima and Kawamura || . It was also 
found there that the critical properties associated with 
the chiral-glass transition were different from those of 
the Ising SG, and that the chiral-glass ordered state ex- 
hibited a one-step- like novel RSB. 

In the chirality scenario of Ref. ||, experimental SG 
transition in real Heisneberg-likc SG magnets is regarded 
essentially as a chiral-glass transition "revealed" via the 
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magnetic anisotropy. Weak but finite random magnetic 
anisotropy inherent to real magnets "recouples" the spin 
to the chirality, and the chiral-glass transition shows up 
as an experimentally observable spin-glass transition. 

The purpose of the present Letter is to reexamine the 
SG ordering of the 3D isotropic Heisenberg SG under ap- 
plied fields in light of the above chirality scenario. We 
first argue some of the possible consequences of the chiral- 
ity picture on the finite-field SG properties, and then per- 
form extensive Monte Carlo (MC) simulations to check 
how the spin and chirality really order in fields. We show 
that the chiral-glass phase, essentially of the same char- 
acter as the zero-field one, remains stable as a true ther- 
modynamic phase in applied fields. It is also found that 
the associated chiral-glass transition line possesses some 
of the character of the GT line of the mean-field model, 
yet its physical origin entirely different. 

Though we expect that our argument holds quite gen- 
erally, we fix here our model Hamiltonian, which is just 
the one used in our MC simulation below. We consider 
the isotropic classical Heisenberg model on a 3D simple 
cubic lattice, 

n = -J2 J ^i-s j -Hj2st, (i) 

<ij> i 

where H is the intensity of the magnetic field applied 
along the z direction. The nearest-neighbor coupling 
is assumed to take the value J or — J with equal prob- 
ability (±J distribution). Local scalar chirality is de- 
fined for three neighboring spins: Here we define it at 
the i-th site and in the /i-th direction (fi = x, y, z) by 

Xin — "i— 8u • Si x S*i+e M ) being a unit lattice vector in 
the /i-th direction. 

Applied fields reduce the global symmetry of the 
Hamiltonian (1) from the zero-field one Z2 x 50(3) to 
Z2 x 5*0(2), where Z2 refers to the spin reflection with 
respect to an arbitrary plane including the field axis, and 
50(2) to the spin rotation around the field axis. Note 
that, even in fields, the chiral Z2 reflection symmetry, 
characterized by the mutually opposite signs of the chi- 
rality, is kept intact. The chirality is a pseudoscalar in- 
variant under 50(2) (or 50(3) in case of H = 0) rota- 
tions but changes sign under Z2 reflections. The same 
chiral Z2 x 50(2) symmetry also appears in the XY SG 
with two-component spins 0. 

Since the chiral Z 2 is supposed to be decoupled from 
the spin 50(3) already in zero field, and the applied field 
serves only to reduce the decoupled 50(3) to 50(2), one 
naturally expects that the Z2 chiral-glass transition, es- 
sentially of the same type as the zero-field one, would 
persist even under magnetic fields fl0|| . More specifically, 
the chiral-glass transition in a field should lie in the same 
universality class as the zero-field one, characterized by 
the same set of exponents, and the chiral-glass ordered 
state in a field exhibits the same kind of one-step-like 
RSB as in zero field. 

Since the chiral-glass order in a field is expected to 



be essentially of the same character as in zero field, the 
chiral-glass transition temperature under magnetic fields 
Tcg(H) should be a regular function of H. If one takes 
account of the obvious symmetry H «-> —H, the chiral- 
glass transition line at low enough fields should behave 
as 

T CG (0) - T CG (H ) = cH 2 + c'H 4 + • • • . (2) 

Generally, the coefficient c could be either positive or 
negative. Interestingly, the above form of the chiral tran- 
sition line is similar to the so-called GT line of the mean- 
field Heisenberg SK model. We emphasize, however, that 
their physical origin is entirely different. The quadratic 
dependence of the chiral-glass transition line is simply 
of regular origin, whereas that of the GT-line in the SK 
model cannot be regarded so. 

Concerning the spin order, applied fields trivially in- 
duce a nonzero longitudinal (parallel to the field) spin 
order at any temperature. The behavior of the trans- 
verse (perpendicular to the field) component could be 
more nontrivial. If one recalls the recent numerical re- 
sults on the 3D XY SG which indicate the absence of 
the standard SG order at least just below the chiral-glass 
transition [p||Tl|l , the transverse spin order is also unlikely 
to arise in the present case, at least just below the Z2 
chiral-glass transition. 

In order to examine whether the above expectation 
really holds or not, we next perform extensive MC sim- 
ulations on the isotropic ±J Heisenberg SG model (1). 
Simulations are performed for a variety of fields H / J = 
0.05, 0.1, 0.5, 2.0, 3.0, 5.0. The lattices studied are simple- 
cubic lattices with L 3 sites with L = 6 ~ 16 with periodic 
boundary conditions. The system is fully equilibrated 
with use of the temperature-exchange method [Q. We 
show here the data explicitly for the particular field value 
of H/J = 0.5, where the sample average is taken over 
128-800 bond realizations. 

By running two independent sequences of systems 
(replica 1 and 2) in parallel, we compute a scalar chi- 
ral overlap q x between the chiralities of the two replicas 

b y 1x = Jn X«m Xi^Xi^i as wel1 as a spin-overlap ten- 
sor q^u between the transverse n and v components (/z, 

v=x,y) of the spin by q^ = jfJ2i S il S i?- Then, in 
terms of these overlaps, we calculate the Binder ratios of 
the chirality g x , and of the transverse {XY) components 
of the spin gx defined in the standard manner: See Ref. 
H for detailed definition. The results are shown in Fig.l. 
The Binder ratio of the chirality g x exhibits a negative 
dip which, with increasing L, tends to deepen and shift 
toward lower temperature. Furthermore, g x of various L 
cross at a temperature slightly above the dip temperature 
T<jip on negative side of g x , eventually merging at tem- 
peratures lower than Tdi p . The observed behavior of g x is 
similar to the one observed in zero field . As argued in 
Ref. H , the persistence of a negative dip and the crossing 
occurring at g x < 0, are strongly suggestive of the occur- 
rence of a finite-temperature transition where 3x(-^cg) 
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and g x (Tcc) take negative values in the L — > oo limit. 
In the inset of Fig. 1(a), we plot the negative-dip tem- 
perature Tdip(-L) versus 1/L. The data lie on a straight 
line fairly well, and its extrapolation to 1/L — gives an 
estimate of the bulk chiral-glass transition temperature, 
Tgg/J ~ 0.25. (More precisely, Tdi P (L) should scale with 
L x l v where v is the chiral-glass correlation-length expo- 
nent. As shown below, our estimate of v ~ 1.3 comes 
close to unity, more or less justifying the linear extrap- 
olation employed here. Extrapolation with respect to 
£i/i-3 yields T CG /J ~ 0.23.) 

In sharp contrast to g x , Binder ratio of the transverse 
component of the spin gT decreases monotonically to- 
ward zero with increasing L, without a negative dip nor 
a crossing, suggesting that the transverse component of 
spin remains disordered even below Tqg- 




FIG. 1. The temperature and size dependence of the 
Binder ratios of the chirality (a), and of the transverse (XY) 
component of the spin (b), in a field H/J = 0.5. Inset dis- 
plays the negative-dip temperature vs. 1/L. The solid and 
broken lines are the best fits assuming the 1/L and 1/L 1 / 1,3 
dependence, respectively. 

Validity of our estiamte of Tgg has been checked also 
from the behavior of other quantities, e.g., the temporal- 
decay of the equilibrium chirality autocorrelation func- 
tion (data not shown here), which yields Tgg I J = 
0.23(2) in accord with the above estimate. In order to 
probe the possible RSB in the chiral-glass ordered state, 
we display in Fig. 2 the distribution function of the chiral- 
overlap defined by P(q x ) = [{5{q x — q' x ))] where < • • • > 



and [• ■ ■] represent the thermal average and the sample 
average, respectively. The existence of a growing "cen- 
tral peak" at q x = for larger L, in addition to the stan- 
dard "side-peaks" corresponding to ±?cG' suggests the 
occurrence of a one-step-like peculiar RSB in the chiral- 
glass ordered state. Similar behavior was observed in 
the chiral-glass state in zero field || . The existence of a 
persistent negative dip in the Binder ratio g x is also con- 
sistent with the occurrence of such a one-step-like RSB 
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FIG. 2. Chiral-overlap distribution function at 

T/J = 0.16, well below the chiral-glass transition point 
Tcg/J ^ 0.23, in a field of H/J = 0.5. 

With setting Tqg/J = 0.23 as determined above, we per- 
form the standard finite-size scaling of the chiral-glass or- 
der parameter [< q x >] and of the chiral autocorrelation 
function, to estimate various chiral-glass exponents. We 
then get v = 1.3(2), r\ = 0.6(3), z — 5.3(5), which turns 
out to agree within errors with the corresonding zero-field 
exponents of Ref. [^J. The results seem consistent with 
a common universality class occurring both in zero-field 
and finite-field chiral-glass transitions. 

Similar calculations and analysis are repeated for other 
field values as well. For H/J = 0.1, in particular, we have 
performed the same scale of intensive calculation as was 
done for H/J = 0.5, to find that all qualitative features 
are similar. The chiral-glass transition with a one-step- 
like RSB occurs at Tqg/J — 0.21(2), with the exponents 
v = 1.3(2), T] — 0.6(3), z — 4.9(5), which agrees within 
errors with our estimates for H/J = 0.5. 

By collecting the Tqg values determined for other field 
values, we construct a phase diagram in the temperature- 
magnetic field plane as shown in Fig. 3. The chiral-glass 
state remains quite robust against magnetic fields. In- 
deed, Tgg(H) is not much reduced from the zero-field 
value even at field as large as ten times of Tcg(0). This 
somewhat surprising property probably arises from the 
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fact that the magnetic field couples in the Hamiltonian 
directly to the spin, not to the chirality, and the effec- 
tive coupling between the field and the chirality is rather 
weak. At lower fields, the chiral-glass transition line is 
almost orthogonal to the H = axis, consistent with 
the expected behavior Eq.(2). Our data are even not in- 
consistent with the coefficient c in Eq.(2) being slightly 
negative so that Tqg{H) initially increases slightly with 
H, though it is difficult to draw a definite conclusion due 
to the scatter of our estimate of Tqg (H) . 
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FIG. 3. The temperature - magnetic field phase diagram 
of the 3D isotropic ± J Heisenberg spin glass. Note the differ- 
ence in energy scales of the temperature- and the field- axes. 



observation of the field-insensitive chiral-glass transition 
line JiJUJlJ. Other interesting possibility revealed by our 
analysis is that the SG transition line might extend to 
higher fields than hitherto suspected. For example, for 
canonical SG AuFe, Campbell et al recently determined 
the SG phase boundary in fields by torque measurements 
0: At lower fields, applied fields rapidly suppress the SG 
transition giving rise to the standard AT-like behavior, 
whereas around the maximal fields of the measurements 
(~ 7T) the SG transition line becomes almost-field inde- 
pendent. If one roughly estimates the effective J of this 
system from its zero-field SG transition temperature, it is 
of order 50K. If one assumes that the energy scale of our 
present model calculation could roughly be applicable to 
AuFe, the SG ordered state should extend to fields much 
higher than I0T without much reduction, or even with 
a slight increase, in T g , although, considering the differ- 
ence in microscopic details between the present model 
and real AuFe, one cannot expect a truely quantitative 
correspondence. Anyway, further high-field experiments 
on AuFe and on other Heiscnberg-like SG magnets might 
be worthwhile to determine the SG phase boundary in 
the high-field regime. 

The numerical calculation was performed on the Hi- 
tachi SR8000 at the supercomputer center, ISSP, Uni- 
versity of Tokyo. The authors are thankful to Dr.K. 
Hukushima for useful discussion. 



Finally, we wish to discuss implications of our re- 
sults to real experimental SG. In real Heisenberg-like SG 
magnets, weak but finite random magnetic anisotropy 
neglected here plays a role. Among other things, the 
anisotropy recouples the spin to the chirality, trans- 
forming the chiral-glass state into the spin-glass state. 
Furthermore, in the presence of both random magnetic 
anisotropy and magnetic field, all global symmetries of 
the Hamiltonian will be lost. Nevertheless, we expect 
that the chiral-glass transition should still persist in fields 
as a pure RSB transition, not accompanying the global 
.^-symmetry breaking. Difference in the broken sym- 
metries in zero- and in finite fields, however, causes a 
singular crossover behavior in fields, in sharp contrast 
to the fully isotropic case. We believe that this singu- 
lar crossover line expected in the weak-field regime of 
anisotropic system is nothing but the AT-like transition 
line ubitiously observed experimentally. Further detailed 
nature of this crossover expected in the low-field regime 
will be discussed in a separate paper Q . 

Meanwhile, in the strong-field regime where the ap- 
plied field dominates the anisotropy, main features of the 
experimental phase diagram would be described by our 
present results on isotropic system. Thus, the experimen- 
tal observation of the GT-like field-insensitive transition 
line in the high-field regime is fully consistent with our 
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